Dressed Counterions: Poly- and Monovalent Ions at Charged Dielectric Interfaces 
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We investigate the ion distribution and overcharging at charged interfaces with dielectric inho- 
mogeneities in the presence of asymmetric electrolytes containing polyvalent and monovalent ions. 
We formulate an effective "dressed counterion" approach by integrating out the monovalent salt de- 
grees of freedom and show that it agrees with results of explicit Monte-Carlo simulations. We then 
apply the dressed counterion approach within the framework of the strong-coupling theory, valid 
for polyvalent ions at low concentrations, which enables an analytical description for salt effects as 
well as dielectric inhomogeneities in the limit of strong Coulomb interactions on a systematic level. 
Limitations and applicability of this theory are examined by comparing the results with simulations. 



I. INTRODUCTION 

Coulomb interactions are important in almost all soft 
matter systems at room temperature Salient fea- 
tures of Coulomb systems such as ionic solutions arise 
as a consequence of the extreme long range of these in- 
teractions. In the last decade, concerted efforts were 
invested in order to understand the statistical behav- 
ior of soft matter systems in the presence of polyva- 
lent ions. It became evident that such ions can be- 
have very differently from monovalent ones 043 ■ The 
reason lies in different strengths of electrostatic interac- 
tions between constituents of multi-component ionic so- 
lutions. The distance at which the interaction energy 
between g-valcnt ions equals the thermal energy k^T is 
g^^B, where £b = eQ/(47reeofcB'r) is known as the Bjer- 
rum length (with cq being the elementary charge and and 
e the dielectric constant of the solvent). Similarly, the 
distance at which each ion interacts with an oppositely 
charged surface (of absolute surface charge density a) 
with an energy scale k^T is the Gouy-Chapman length, 
/.t = e^l {iirqt-Qa). The behavior of a charged system 
depends on the ratio of the two, which defines an electro- 
static coupling parameter S ~ q^l^/ ^ @. In many prac- 
tical cases for monovalent ions this coupling parameter 
is 5 < 1, which justifies a mean-field Poisson-Boltzmann 
(PB) description of the system. In contrast, for polyva- 
lent ions 5^1, which leads to diverse phenomena, such 
as overcharging and like-charge attraction ■ 

A strong-coupling (SC) theory valid for S ^ 1 pro- 
vides a consistent and systematic theoretical framework 
for understanding Coulomb systems in this limit. Until 
now it has been formulated explicitly for a counterion- 
only system Q, which has been extended to a case with 
dielectric inhomogeneities as well [l2|. We now make 
a crucial step forward (in order to study the more re- 
alistic and experimentally relevant cases) and formulate 
the SC theory for a mixture of symmetric and asymmet- 
ric electrolytes (poly- as well as monovalent salts) in the 



presence of dielectric inhomogeneities. 

In order to investigate the salient features and to test 
our theoretical approach to highly asymmetric multicom- 
ponent Coulomb systems we set ourselves to study the 
polyion concentration profiles and the possibility of over- 
charging by three widely different approaches. The first 
one is the standard (explicit) Monte Carlo (MC) simu- 
lation, which treats all ions explicitly, and will merely 
serve as a test for the other two approaches. In the 
second, dressed counterion approach, we integrate out 
the monovalent ion degrees of freedom by replacing the 
Coulomb interaction between remaining charges with an 
effective screened Deby e-Huckel (DH) interaction (as first 
described in Ref. [Tj] for a dielectrically homogeneous 
system with uniform screening). Such screened interac- 
tions arc then used to perform implicit MC simulations 
that treat all charges explicitly except the monovalent 
ions, and show an extremely good agreement with ex- 
plicit MC, confirming the validity of the approach based 
on dressed counterions in highly asymmetric systems 
q ^ 1. Finally, we im plem ent the dressed counterions 
strong coupling theory [l3|, which allows us to obtain 
an analytical description of ionic distributions and over- 
charging in the limit of strong coupling with the local 
surface field. The validity of the SC dressed counterion 
theory is thus examined by making a comparison with 
both explicit and implicit simulations. 



II. THE MODEL 

The model system considered here. Fig. [U comprises 
an impenetrable charged planar dielectric interface with 
surface charge density —a bathed in a mixture of mono- 
valent 1:1 and polyvalent q:\ salts with bulk concentra- 
tions of no and cq, respectively. The total bulk concen- 
tration of monovalent (-1-1) counterions is therefore no, 
that of polyvalent {+q) counterion concentration Co, and 
monovalent (—1) colon concentration no +gco. We shall 
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FIG. 1: (Color online) Schematic view of an impenetra- 
ble charged dielectric interface of surface charge density —a 
bathed by a monovalent 1:1 salt and polyvalent q:l salt, with 
larger q-valent counterions. For water-hydrocarbon interface 
one has e = 80, e' = 2, which gives dielectric A = 0.95 that 
allows for an approximate treating as A = 1. 



refer to +q polyvalent counterions simply as "counteri- 
ons" , which should not be confused with monovalent salt 
counterions. We will express the bulk concentration of 
polyvalent counterions by the parameter 



X 



(1) 



The dielectric constant of the surface e' can be sig- 
nificantly different from the solvent e. Wc define the 
dielectric contrast parameter as 



A = 



e — £ 



(2) 



The presence of a large dielectric contrast at the interface 
(A ^ 1, or e ^ e') can lead to enhanced effects due to 
the presence of polyvalent counterions @, [H, H, Eli • 

Note that no overcharging would be possible if the solu- 
tion consisted only of neutralizing counterions. However, 
this situation is expected to be modified significantly by 
the presence of monovalent salt, where various degrees of 
overchargingof the charged surface have been reported 
previously |16l - [2l| . 



III. DRESSED COUNTERION THEORY 

In what follows we will apply the "dressed counterion" 
approach to our model system, which was introduced 
recently [Tsj in a different context and was not thor- 
oughly examined by explicit simulations. This approach 
is based on the fact that for highly asymmetric solutions 
with q ^ 1, the monovalent ions can be integrated out 
from the partition function and their effect is replaced 
by an effective ("dressed") DH interaction that acts be- 
tween the remaining polyions and the surface charges. 
The idea stems from a realization that monovalent ions, 
because of their small valency, are only weakly coupled 
to all other charges, i.e., polyvalent counterions as well 
as surface charges. As a result, they can be integrated 



out, a procedure that can be justified systematically in a 
full field-theoretical framework [l^ ■ 

For the sake of simplicity, we shall illustrate this idea 
in the context of mean-field theory where the system is 
described via the Poisson-Boltzmann equation. 



'AniB [n-oe ^' - (no -I- qco)^^ + qc^e' 



(3) 

where the first three terms correspond to -1-1 ions, —1 
ions and +q ions, respectively. The fourth term /5o(r) 
stands for the external surface charges (in units of eg) 
and V' is the dimensionless mean electrostatic potential 
(in units of fcer/ep). In the situation when -0 1 and 
q ^ 1, one can linearize the first two terms in the above 
PB equation as 

VV = - 47r^B [gcoe-"^ - gco + po(r)] . (4) 

The constant zcroth-ordcr term, — gcg, can be absorbed 
into the definition of the electrostatic potential ij:, corre- 
sponding to a displacement of the potential by a constant 
Donnan potential. This procedure changes nothing for 
our subsequent analysis. We have hence introduced the 
screening parameter 



(5) 



which is defined by the total amount of all monovalent 
ions, i.e., ions coming from the 1:1 salt with concen- 
tration no as well as the additional monovalent colons 
that enter the solution via the polyvalent q:l salt with 
concentration qcQ. This form of the effective screening 
parameter follows also from more general statistical me- 
chanical grounds in the context of colloidal interactions 
in the presence of low concentration of salts [23- It is 
also amended compared to the definition in our previ- 
ous work in Ref. [l3| , where the regime of relatively high 
monovalent salt concentration was of main interest. 

Note however that the above PB equation cannot prop- 
erly treat polyvalent ions due to high correlation effects. 
We now make a crucial step, namely we only treat mono- 
valent ions in the linearized PB sense, whereas polyvalent 
ions along with surface charge are considered explicitly, 
expressed in terms of explicit charge density poxp(r), viz. 



47r£BPoxp(i")- 



(6) 



This linearized PB equation then leads straightforwardly 
to the relevant Green's function or the non-homogeneous 
DH kernel, representing the screened DH interaction 
(modified by the presence of any dielectric interfaces) 
between all the remaining explicit charges after the sim- 
ple salt degrees of freedom are integrated out. It can 
be represented for the present system via a 2D Fourier 
transform as 



uy)b{Q\z,z') 
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ek — e'' 
ek + e'' 



(7) 
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where Q is the transverse wave-vector and k = 
a/k^ + Q^. The first term is the usual (homogeneous) 
DH interaction in free space, and the second term is the 
non-homogeneous correction due to the presence of di- 
electric discontinuities and the inhomogeneous distribu- 
tion of salt ions in the system (leading to screening dis- 
continuities, i.e., no salt ions are present in the region 
z < -see Fig.[T|). 

Now in order to calculate the full partition of the sys- 
tem in any given setup (beyond the mere mean-field limit 
approximation that we used for the illustration above), 
one has to integrate over the degrees of freedom associ- 
ated with the explicit polyvalent counterion charges as 
well (see Ref. [l^ for a systematic field-theoretical treat- 
ment). We shall consider this procedure later for the 
study of the system in the strong-coupling limit. 

The dressed counterion approach thus enables an enor- 
mous reduction in the number of explicit particles in 
the system, comprised now only of polyvalent counteri- 
ons and fixed surface charges that interact via "dressed" 
non-homogeneous DH interactions. Although substan- 
tially simplified, the "dressed counterion" model still rep- 
resents a full many-body problem (that can be used to 
study the system in all regimes of the coupling parameter 
by prescribing the modified pair interaction according to 
Eqs. ([5]) and ©) and is as such exactly solvable only by 
computer simulations. 

In many biologically relevant cases the dielectric ma- 
terial of the bounding surfaces has a dielectric constant 
e' ~ 2 — 4, as is the case with e.g. hydrocarbons or 
DNA hydrophobic core, whereas the solvent is water with 
e ~ 80. The notion that e ^ e' (A ~ 1) can sim- 
plify theoretical treatment enormously, since the non- 
homogeneous DH kernel, Eq. ([7]), can in that case be 
rewritten as 



UDYi{Q\Z,z') 
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(8) 



or in real space simply as 



UDHlr, r') 



-K |r— r' I 



-K,\r—v'--\-2z'Gz I 
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(9) 



This brings us to the concept of screened image charges, 
where the presence of each charge produces its own im- 
age on the other side of the interface (second term). Here 
Bz = (0, 0, 1) stands for unit vector perpendicular to the 
surface. In this limit the exact value of e' does not en- 
ter explicitly. If not stated otherwise we will mostly deal 
with the case where A ~ 1, which is also a very good ap- 
proximation to the realistic A = 0.95 case as will be seen 
later. Fig. [3]:. Note that the above form of the Green's 
function takes into account the standard dielectric im- 
ages as well as the ionic cloud images stemming from the 
inhomogeneous distribution of salt in the system. 



A. Implicit vs. explicit simulations 

To examine the validity of the dressed counterion ap- 
proach (for a highly asymmetric system with a large q) 
we first directly compare the results of implicit and ex- 
plicit MC simulations. In order to perform the MC sim- 
ulations with the dressed counterion ansatz that we also 
refer to as the implicit MC simulations, we need to evalu- 
ate the interactions between all the charged via the non- 
homogeneous DH kernel, Eqs. ([8]) and (|9]). 

The interaction between every dressed counterion and 
the surface charge is then written as 



Oc 



13 J po{r)uj:,H{r,r')pc{r')drdr' 



(10) 



with surface charge operator pair) ~ —a5{z) and counte- 
rion charge operator Pci^) = eo(7(5(R— r) with R standing 
for the counterion position. 

The interaction between dressed counterions is addi- 
tive. For every pair of different counterions i and j, 
(i ^ j), the interaction is 



PWcc{i,j) = /3(eo(7)^itDH(r»,rj). 



(11) 



Each counterion itself also experiences a self-energy due 
to the dielectric half-space. 



4z 



(12) 



This repulsive interaction is responsible for the depletion 
effect in the vicinity of the surface. Its range is twice as 
short as attractive surface-charge interaction, Eq. (|10p . 
The dielectric self-energy of the ion represents a general- 
ization of the concept of Born energy valid for an ion in a 
homogeneous dielectric. In the present paper we do not 
take into account the finite size of the ions in the dielec- 
tric self-energy, but we do take into account the presence 
of the dielectric inhomogeneities. 

At this stage we introduce the dimensionless (rescaled) 
units by defining 



r/p, k = pK, X = PX- 



(13) 



In order to connect with unsealed values, table |T] gives 
illustrative examples of actual concentrations for vari- 
ous dimensionless parameters k and x for S = 50 and 
S = 100. Although higher values of dimensionless screen- 
ing K > 0.5 can be difficult to achieve in aqueous solutions 
due to the high salinity required, it can be nevertheless 
more easily achieved in solvents with lower dielectric con- 
stant (e.g., ethanol). 

In our second approach we treat all the ionic species 
on the same footing by explicitly including them into 
MC simulations which we then refer to as the explicit 
MC simulations. For demonstration purposes we fo- 
cus on the (7 = 4 counterions and a surface charge 
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H = 50 

K = 0.1 0.2 0.5 0.8 


H = 100 
K = 0.1 0.2 0.5 


X = 0.1 
0.2 
0.5 
1.0 


no = 15mM 68mM 440mM I.IM 
9mM 62mM 430mM I.IM 
- 15 mM 390 mM I.IM 
- 220 mM 0.9 M 


Co = 1.1 mM 
4.7mM 
30 mM 
120 mM 


no = 63 mM 280 mM 1.8 M 
36 mM 250 mM 1.8 M 
- 63 mM 1.6 M 
- 0.9M 


Co = 4.5 mM 
18 mM 
UOmM 
450 mM 



TABLE L Actual values of concentrations no and Co in aqueous environment as follow from rescaled parameters k and x for 
tetravalent counterions {q = 4) and surface charge densities a = 0.25eo/nm^ (H — 50) and a = 0.5eo/nm^ (H = 100). Note 
that Co depends only on Xj i-e-, co = co(x), whereas no = no(x, fi) depends both on x and k. 



density cr = 0.25eo/nni^. In this case fi — 0.23 nm 
whereas S = 50. In contrast, the coupling parameter 
for monovalent ions is Smono = 0.8. To avoid Bjer- 
rum pairing [23| we have used closest-approach distances 
~ 0.35 nm between monovalent cations and an- 
ions, and (i+4^_i = 0.45 nm between (polyvalent) counte- 
rions and monovalent colons. 



IV. RESULTS FROM IMPLICIT AND 
EXPLICIT SIMULATIONS 

The results for (polyvalent) countcrion density profiles 
obtained from both types of MC simulations are shown 
in Fig. [5^1- Counterions arc localized near the charged in- 
terface but arc depleted from the immediate vicinity due 
to image repulsions. At large distances the profile satu- 
rates to the bulk value, Cq. The higher the monovalent 
salt concentration, i.e., the larger the screening k, the 
smaller the amount of counterions close to the surface as 
larger amounts of salt screen the surface charge more ef- 
ficiently. The comparison of explicit (dots) and implicit 
(solid line) MC simulations shows not only qualitative 
but also good quantitative agreement. 

An interesting quantity to examine is the cumulative 
charge, Q{z), being the sum of all charges within a finite 
distance from the surface, including the surface charges 
themselves 

1 r 

Q{z) = - [eon+{z') - ean^iz') 
^ Jo 

+ eoqciz')-poiz')]dz', (14) 

where n± (z) represent the density of monovalent ions and 
c{z) that of polyvalent counterions. Note that the Q{z) 
is normalized so that it reaches Q = —1 at the surface 
{z = 0) and approaches far away from the surface (z — >■ 
oo) due to complete neutralization of the surface charge. 

In the case where the ions are treated explicitly the 
evaluation of Q{z) is straightforward. Within the im- 
plicit dressed counterion approach the monovalent salt 
densities in the definition of Q need to be recalculated 
from the DH equation, viz. n+(r) — n_(r) ~ —2nQtp{r), 
where ijj corresponds to the dimcnsionlcss electrostatic 
potential that is evaluated from known explicit densities. 



VIZ. 



V'(r) = /3eo j UDH(r, r')[eogc(r') + po(r')] dr'. (15) 

Inserting the surface charge density po{v) = —aS{z) and 
the non- homogeneous DH kernel, Eq. ([5]), we get after 
some algebraic manipulation the expression for the cu- 
mulative charge in dimensionless units 



Qiz) - -c-^^ + -x^ 



(16) 



{z - z')e-«l^-^'l + e-'=(^+^')l d{2') dz' 



sgn 



The first term is due to the DH exponential screening 
in the absence of polyions, whereas the second one cor- 
responds to the effects due to polyvalent ions, where 
c(z) = c(z)/co. 

As seen in Fig. [2]3 the surface charge is neutralized 
more quickly in the case with a higher salt screening. 
If the screening is low enough, the surface can be over- 
charged leading to positive values of Q at a finite z. The 
overcharging degree turns out to be relatively small (in 
fact an order of magnitude smaller) compared to the mag- 
nitude of the bare surface charge, which agrees quantita- 
tively with results from several recent explicit MC studies 
on similar systems [iB, [13, HH but disagrees with oth- 
ers [Mil- 

The agreement between implicit and explicit models is 
less pronounced close to the surface as the DH approach 
used in implicit simulations does not account for finite 
size of monovalent ions as well as image effects in the dis- 
tribution of monovalent ions. Such discrepancies of ionic 
distributions in close vicinity to the interface however 
do not affect the cumulative charge at larger distances, 
where implicit and explicit results match perfectly. 

We shall quantify the discrepancy between implicit 
and explicit MC results by comparing the relative dif- 
ference in the surface areas under the counterion profile 



curves. We thus calculate the integrals AS = J \ 



-imp 



(^)- 



Ccxpiz)\dz and S = J Cexp{z)dz within the interval corre- 
sponding to 2xFWHM (full width at half maximum) on 
both sides of the peak. The relative error is then defined 
as err = AS/S, and is presented in Fig. [21:;. Note that ac- 
cording to the definition of k we should have x ^ 
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FIG. 2: (Color online) Comparison of explicit and implicit MC simulations for S = 50. (a) Counterion density profile at fixed 
bulk amount of counterions x and two different salt concentrations k. Implicit MC results (solid line) compare quite well with 
explicit MC results (dots), (b) Cumulative charge that follows from the results in (a). In the case of smaller screening {k — 0.30) 
the cumulative charge reaches positive values, which corresponds to overcharging of 3%. (c) "Phase diagram" showing the 
deviations between implicit and explicit MC simulations for H = 50. The agreement is better for smaller values of k. 



As can be seen in Fig. ^ the agreement between implicit 
and explicit MC results is fairly good in a wide range of 
screening parameters. The agreement becomes less sat- 
isfying at large screenings and to some extent also when 
X is increased, i.e., at larger concentrations of the poly- 
valent counterions. 



These results may seem in contradiction with a naive 
expectation that the DH approach performs better for 
larger screenings. Note however, that the finite size of 
monovalent ions becomes increasingly more important 
at higher screenings (salt concentrations), an effect not 
taken into account within the present DH approach. On 
the other hand, one might also expect a complete failure 
of the DH approach at low screenings, i.e., when k < 
especially near the charged surface where non-linear ef- 
fects become important. In this case it turns out that the 
dielectric mismatch acts in the opposite direction and the 
image repulsion depletes the counterions from the imme- 
diate vicinity of the surface, hence the exact form of the 
interaction with the surface charge in its vicinity is not 
important. Yet at very low screenings, some additional 
discrepancies might be observable due to non-linear ef- 
fects. However, down to k ~ 0.1 we still find very good 
agreement within err < 10%, c.f. Fig. 



From Fig. [2t it thus appears that the dressed coun- 
terion theory is a good approximation for monovalent 
salt mediated interactions. This remains true in a very 
large segment of the parameter space where the discrep- 
ancies seldom exceed 20%, clearly corroborating the prac- 
tical value of the dressed counterion approach (for highly 
asymmetric cases with g 3> 1) and the implicit simula- 
tions, which can be orders of magnitude less expensive 
than the standard explicit MC simulations, especially for 
very large systems of monovalent and polyvalent salt mix- 
tures. 



V. SC DRESSED COUNTERION LIMIT 

We now focus on the analytical SC description of our 
system. As already stated, the monovalent and polyva- 
lent ions need separate treatments due to their vastly 
different coupling to external charges. In the case of 
counterion-only systems the polyvalent counterions can 
be treated in the SC limit S ^ 1 via a virial expansion in 
terms of their fugacity A from the grand-canonical parti- 
tion function Zq as 

zg = 4°' + ^4'^ + ■ ■ ■ ■ (17) 

Here the first term corresponds to the interaction of bare 
surface charges in the absence of any counterions, and the 
second term corresponds to the contribution of a single 
counterion in the system, etc. This approach turns out 
to be asymptotically exact for high values of the coupling 
parameter S . In its original version the theory does 
not include any monovalent salt ions, being limited only 
to the effects of neutralizing counterions at charged sur- 
faces [1]. Here we implement the same idea but take it 
a step further with "dressed" counterions to account for 
asymmetric systems containing additional salt (as in our 
previous work where two charged surfaces were studied 
without any dielectric mismatches p^). One can then 
follow a similar virial expansion approach as above and 
it turns out that the only change in the SC formulation 
wrought by the dressed counterions is to replace the stan- 
dard Coulomb kernel in the original formulation by 
the non-homogeneous DH kernel, Ec^. ([5]), in the dressed 
counterion formulation. This is then the essence of the 
SC dressed counterion theory. Note that in our case we 
focus on a grand-canonical system in equilibrium with 
a bulk solution (rather than a canonical system consid- 
ered in the absence of salt screening 0411)' ^'^d hence 
the fugacity A is explicitly determined by the bulk con- 
centration of (polyvalent) counterions cq. 

From the first-order grand-canonical expansion, 
Eq. P?|) . the grand potential /3<I> = —InZc of the dressed 
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counterions can be evaluated as 



/3$ 



r(0)- 



(18) 



As already mentioned the first-order term stands for 
the one-particle dressed counterion contribution, hence 



Z(j°^ / exp(-/3VFoe(r) - /3W,cif(r)) dr. (19) 



The average number of dressed counterions is then ob- 
tained from the thermodynamic relation 



N = -A(9/3$/9A). 



(20) 



in complete analogy with the standard SC formulation 
for "bare" counterions. The dressed counterion density 
profile can be extracted from the form of the first-order 
term, Eq. 1191 and can be written explicitly 



c{z) 



Co exp — e 



45 



(21) 



again analogous to the expression for the bare counteri- 
ons. The first term on the r.h.s. corresponds to attrac- 
tive surface-counterion interaction, Eq. PH)) . and decays 
exponentially with the screening parameter k. The sec- 
ond term is due to repulsive self-energies of counterions, 
Eq. (|12p. having their origin in counterion self-image in- 
teractions. It decays twice as fast as the first one and 
becomes negligible at larger separations compared with 
the first term. However, at small separations it diverges 
as 1/z and hence engenders a depletion layer of coun- 
terions. Its strength increases with increasing coupling 
parameter S as is known from other image-interaction 
studies [Hill. 

At this stage we should emphasize that our SC dressed 
counterion limit, which by definition corresponds to the 
leading order in the virial expansion, does not necessarily 
correspond to a laterally correlated counterion cloud as is 
the case for the standard counterion-only systems. This 
fact should be taken into account when the validity of 
the SC limit will be considered, see below. 

In Fig. Owe show the main SC results and their trends 
in terms of all three adjustable parameters, 5, k, and A. 
As seen from Fig. increasing the coupling parameter 
S results in a large image repulsion, Eq. (fT2)) . and hence 
a lower density at the peak. Similarly, the density peak 
decreases also by increasing the screening parameter k 
as shown in Fig. ISb. Thus the larger the screening, the 
weaker would be the interaction between all the explicit 
charges and therefore a smaller amount of counterions is 
attracted towards the surface. 

Finally, we also show the density profiles for various 
values of the dielectric contrast A, which requires an ex- 
act treatment of the non- homogeneous DH kernel, Eq. ([7]) 
(rather than the simplified form ([H]))- The counterion- 
surface term Wqc, Eq. (jlO[) . has the same form also when 



the exact kernel is taken into account, and is hence inde- 
pendent of the dielectric contrast. In contrast, the self- 
energy interaction Wgeit takes a more complicated form 
for arbitrary A, viz. 



/3W,eif(r) 



ek - e'Q QdQ 

— — I 

ek + e'Q k 



-2kz 



(22) 



with k = yk'^ + Q'^. As can be seen from numerical 
results in Fig. [3J;, increasing the dielectric contrast re- 
sults in an increased depletion and peak lowering, sim- 
ilar to the case where S or k are increased. Note also 
that in the case of a vanishing dielectric contrast, A = 
or e' = e, there remains slight depletion of counterions 
from the surface vicinity, which is due to solvation ef- 
fects of implicit salt. The latter situation corresponds to 
pure solvent on the left side (with k' = 0) and electrolyte 
solution on the right (with k > 0), which causes a free en- 
ergy loss when every solvated counterion approaches the 
interface. It should be noted also that the results with 
the dielectric contrast A = 0.95 [corresponding e = 80 
and e' = 2] almost match with the limiting case of A = 1, 
which justifies our previous approximations of the non- 
homogeneous DH kernel, Eq. ([5]). 

It is instructive to look at the minimal value of the 
counterion concentration c™'" neccessary for overcharg- 
ing. Fig.|3h' shows the minimal x as a function of k (solid 
lines) which is necessary for overcharging. Above these 
lines the surface is overcharged whereas below it is not. 
Note the nearly linear dependence on S with a very mild 
slope for K < 0.6. 



A. SC criterion 

As already mentioned the SC dressed counterion 
approximation, Eq. (|21[) . differs from the standard 
counterion-only SC theory in two important as- 

pects: the short-range nature of the dressed DH interac- 
tion and the grand-canonical ensemble used for counteri- 
ons. These differences engender several important modi- 
fications in the nature and interpretation of results. First 
of all it follows that on increase of the coupling parameter 
S one does not necessary approach the strong coupling 
limit as understood for counterion-only systems. Instead, 
being a single particle theory with variable number of 
particles, it is expected to be valid when the interactions 
between counterions are negligible compared to thermal 
energy. This leads to a validity criterion which can be de- 
termined as follows. The effective interaction energy of 
a counterion at position r with all the other counterions 
can be estimated in a mean-field fashion as 



/3I^-'(r) ~ /3(eoQ)2 / uuii{r,r')c{r')dr' 



(23) 



The highest contribution is expected for r around the 
profile peak, therefore in order to estimate the maximal 



6 



c(z), 



/ 


' \ A= 1 
\ K = 0.5 


— H = 25 
-- H = 50 
■ ■ ■ H = 100 

E = 200 








/ 

/ ' 
/ / 

/ / 
/ / 


/ 






(a) 



(b) 



(c) 



FIG. 3: (Color online) SC dressed counterion profiles for (a) various coupling parameters, (b) various screening parameters, 
and (c) various dielectric contrast parameters. 
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FIG. 4: (Color online) (a) "Phase diagram" showing the minimal amount of polyvalent salt x that is required to achieve 
overcharging (solid lines) within the SC dressed counterion theory. Dotted lines correspond to the situation where fiW^ = 1, 



Eq. (|23p . with same colors correponding to the same parameters as shown on the graph. Parameter space below dotted lines 
should be well described by the SC dressed counterion approach, (b) Rescaled density profile for dressed counterions for H = 50 
and screening k = 0.4. The SC prediction is independent of x sjad gives good results for low enough x- With increasing x 
the implicit MC results (symbols) show gradual deviations from the theoretical prediction, (c) Cumulative charge, Eq. (|16|) . as 
follows from the profiles in (b) for the SC theory (solid lines) and implicit MC simulations (symbols). 



pW^ we choose r such that it correponds to the posi- 
tion of the maximum of c(r). The SC approach, Eq. (PT|). 
is thus expected to be valid when /3W™^ ^ 1. Since 
this interaction grows with cq, the SC approach is sat- 
isfied when X is small enough. We plot x a-s a function 
of k that corresponds to (3W™^ = 1 in Fig. (dotted 
lines). In the regions below these lines the SC theory 
is expected to work well. Here equal colors of solid and 
dotted lines correspond to the same S. Note that solid 
and dotted lines roughly coincide. This brings us to the 
rule-of-thumb conclusion that the SC theory ((2T|) is valid 
below the overcharging regime. With the emergence of 
overcharging the SC theory then starts to break down. 

The above argument can be confirmed by comparing 
the SC prediction with MC simulations, Figs. HJd, St. 
Note that the rescaled density c(z)/co is independent of x 
in the SC theory, Eq. ((2T|) . As expected the comparison 
with the implicit MC data shows growing disagreement 
with the SC dressed counterion theory as x is increased. 
The disagreement stems from the fact that the theory 



overestimates the counterion density as it neglects the re- 
pulsive counterion-counterion interactions and therefore 
predicts a larger number of counterions near the surface 
as compared with simulations. However, the profile con- 
tinues to display the same qualitative features, i.e., a peak 
with a depletion region, for a very wide range of x. The 
corresponding cumulative charge Q{z), shown in Fig.|4}D, 
shows perfect agreement between the SC theory and sim- 
ulations for the cases below the overcharging regime. The 
SC theory starts to break down quantitatively when the 
system approaches the overcharging point. The devia- 
tions in Fig. |3t become significant for x ~ 0.95. This 
general conclusion can be deduced not only from the cu- 
mulative charge but also from the counterion density pro- 
files and the " phase diagram" , which confirms our rule 
of thumb that the SC prediction is valid below the over- 
charging regime, see below. 

Finally we compare the overcharging threshold value 
of X from simulations and from SC dressed counterion 
theory. The diagram in Fig. [S] shows the demarcation 
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FIG. 5: (Color online) "Phase diagram" characterizing the 
overcharging of a charged surface. The threshold values of 
minimal x neccessary for overcharging are obtained by all 
three different methods, i.e., explicit and implicit MC simu- 
lations as well as the SC dressed counterion theory. 



line between the overcharging regime and the absence 
thereof (for S = 50) as follows from all the three differ- 
ent approaches we have used in this study. The implicit 
and explicit MC simulations agree very well for screening 
parameters k < 0.6, which is in accordance with Fig. [5};. 
On the other hand, the results from the SC theory show 
a shift by a constant value towards a smaller threshold 
value for x- The reason is again that the SC theory over- 
estimates the counterion density near the surface, and 
hence predicts a smaller x needed for achieving the over- 
charging. Nevertheless, all three methods predict a clear 
and consistent trend: a nearly linear dependence of x on 
K at the demarcation line. This fact was also observed 
experimentally, i.e., the more salt is present in the sys- 
tem, the higher concentration of polyvalents is needed in 
order to observe overcharging |24| . As already seen in 
Fig. 13^ the SC theory gives a trend that is only weakly 
dependent on S, which appears to indicate an approxi- 
mate proportionality relation as ^ k. In real units 
this translates into a heuristic relation as 



c'r - no/q', (24) 



consistent with the experimentally observed fact that 
countcrions with higher valency q arc more efficient at 
overcharging, which has also been observed experimen- 
tally [2^. Nevertheless one needs to keep in mind that 
ion-specific effects [2^ also make their mark in strongly 
coupled systems. The SC theory based on purely elec- 
trostatic grounds appears to be unable to account for the 
experimentally observed variation in the strength of the 
SC effects for equally charged polyvalent counterions [l^j . 
The ion specific effects in the SC theory have remained 
largely unexplored. 



VI. CONCLUSION 

In our work we present detailed arguments based on 
extensive simulation for the validity of the dressed coun- 
terion approach and its SC limit, based on a model 
for an ionic solution mixture composed of a monovalent 
(1:1) salt and an asymmetric (q'-l) salt in proximity to a 
charged dielectric interface. 

In the first part we have described the general idea of 
dressed counterions, i.e., the coarse-grained electrostatic 
treatment that replaces explicit monovalent- ion degrees 
of freedom by the effective Debye-Hiickel interaction be- 
tween the remaining charges comprised of the polyvalent 
countcrions and the surface charges. This idea is valid in 
highly asymmetric systems with q ^ 1 and stems from 
the fact that monovalent ions are weakly coupled with 
all the other charges. We compared the MC results of 
dressed counterions (referred to as implicit simulations) 
with the standard explicit MC simulations (referred to 
as explicit simulations) where all the ions are treated ex- 
plicitly. Comparing the results for the density profiles 
of polyvalent countcrions we found excellent agreement 
between both types of simulations. 

Generally, counterions are attracted toward the sur- 
face because of direct electrostatic interactions but are 
expelled from the vicinity of the surface due to the di- 
electric image charges and in general solvation repulsion. 
For large polyvalent salt concentration we observed over- 
charging, i.e., the net excess of opposite charge near the 
surface. Both types of simulations agree fairly well on 
the degree of overcharging, which is typically of the or- 
der of few percent of the surface charge, in accordance 
with recent MC studies [H, [13, [HI . 

In the second part of the paper we implemented the 
analytical SC dressed counterion approach, based on a 
strong-coupling one-particle description , to the model 
at hand. The SC dressed counterion theory gives a very 
simple analytical expression for the counterion density 
profile, Eq. (|2ip . and also enables simple predictions for 
the onset of overcharging. We estimated the regime of 
validity of the SC theory based on the fact that the 
counterion-counterion interactions should be negligible 
when compared to the thermal energy. As it turns out 
this criterion approximately coincides with the onset of 
overcharging. In other words, the SC dressed counterion 
theory is valid in the regime just below the overcharg- 
ing regime, which was confirmed by implicit MC simu- 
lations. In its regime of validity the SC dressed counte- 
rion theory gives good predictions for the density profile 
of counterions and can also approximately predict the 
onset of overcharging. In agreement with the MC sim- 
ulations, the theory shows a linear relationsip between 
the polyvalent salt concentration needed for overcharg- 
ing and the monovalent salt concentration in the system. 
It also predicts higher-valency counterions are more ef- 
fective in overcharging the surface. These results agree 
with recent experimental findings [l^, • 

Finally we should emphasize that the merit of the 
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dressed counterion approach is firstly in its practical use 
in implicit MC simulations, which give satisfactory re- 
sults and can be several orders of magnitude less ex- 
pensive compared to explicit MC simulations. Secondly, 
it can be easily included in analytical frameworks for 
charged systems (e.g., the strong coupling theory in our 
case), which may thus allow for a deeper understanding 
of the general features of the charged systems that con- 
tain highly asymmetric electrolytes. 
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